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Abstract 

We consider the elastic energy of a hanging drape - a thin elastic sheet, pulled down by 
the force of gravity, with fine-scale folding at the top that achieves approximately uniform 
confinement. This example of energy-driven pattern formation in a thin elastic sheet is of 
particular interest because the length scale of folding varies with height. We focus on how the 
minimum elastic energy depends on the physical parameters. As the sheet thickness vanishes, 
the limiting energy is due to the gravitational force and is relatively easy to understand. 
Our main accomplishment is to identify the “scaling law” of the correction due to positive 
thickness. We do this by (i) proving an upper bound, by considering the energies of several 
constructions and taking the best; (ii) proving an ansatz-free lower bound, which agrees 
with the upper bound up to a parameter-independent prefactor. The coarsening of folds in 
hanging drapes has also been considered in the recent physics literature, using a self-similar 
construction whose basic cell has been called a “wrinklon.” Our results complement and 
extend that work, by showing that self-similar coarsening achieves the optimal scaling law 
in a certain parameter regime, and by showing that other constructions (involving lateral 
spreading of the sheet) do better in other regions of parameter space. Our analysis uses a 
geometrically linear Foppl-von Karman model for the elastic energy, and is restricted to the 
case when Poisson’s ratio is zero. 


1 Introduction 

We consider a hanging drape - a thin elastic sheet, pulled down by the force of gravity, with fine- 
scale folding at the top that achieves approximately uniform confinement. We are interested in 
how the shape of the sheet varies with height. Since bending costs elastic energy, one expects to 
see less bending far from the top. This effect can be achieved by two rather different mechanisms: 
coarsening of the folds, or spreading of the sides. Our analysis includes a study of these two 
mechanisms and how they interact. 

Our viewpoint is variational: we focus on how the minimum elastic energy depends on 
the physical parameters. Our elastic energy functional has three terms: the membrane energy 
(which penalizes stretching), the bending energy (which penalizes bending), and a loading term 
(representing the effect of gravity). The membrane energy is nonconvex, while the bending 
energy acts as a regularizing singular perturbation. To capture the essential behavior in the 
simplest possible setting, we use a Foppl-von Karman model with Poisson’s ratio 0 for the 
membrane energy. The limiting behavior as the sheet thickness h 0 is then quite easy to 

‘Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany. Email: bella@mis.mpg.de. This 
work was begun while PB was a PhD student at the Courant Institute of Mathematical Sciences. Support from 
NSF grant DMS-0807347 is gratefully acknowledged. 

'Courant Institute of Mathematical Sciences, New York University, New York, NY. Email: 
kohn@cims.nyu.edu. Support is gratefully acknowledged from NSF grants DMS-0807347, DMS-1311833 and 
OISE-0967140. 


1 



understand: the sheet hangs straight down, under tension due to the force of gravity. The 
energy of this deformation is easy to find. Our main accomplishment is to assess (at the level 
of scaling) the leading-order correction due to positive h. We do this by (i) proving an upper 
bound on the minimum energy, by considering several constructions and taking the best; and 
(ii) proving an ansatz-free lower bound, in which the correction due to positive h agrees - up to a 
parameter-independent prefactor - with the upper bound. (This summary is slightly misleading: 
in fact the thickness enters our bounds through a nondimensional ratio that also involves the 
gravitational force, the curtain’s height, and the length scale of the wrinkling imposed at the 
top; moreover our results are not asymptotic as h —>• 0 - rather, they require mainly that h be 
small compared to the length scale of the wrinkling imposed at the top. For a precise statement 
of our bounds see Section [3]) 

The coarsening of folds or wrinkles in a sheet under tension has also been considered in 
the physics literature is iia nasals]. Our work is especially strongly connected with the 
treatment of “heavy sheets” in m which uses a self-similar construction whose basic cell is 
called a “wrinklon.” One of the constructions used for our upper bound involves coarsening 
of the folds with no lateral spreading; it is essentially the same as the deformation studied in 
m (see Remark EH). Our results complement and extend those of that paper, by showing 
that the self-similar coarsening considered there achieves the optimal scaling law in a certain 
parameter regime, and by showing that other constructions (involving lateral spreading of the 
sheet) do better in other regions of parameter space. We note in passing that |3T| also discusses 
the coarsening of folds in “light sheets.” An analysis of that problem in the spirit of the present 
paper has (to the best of our knowledge) not yet been done, though the blistering problem 
considered in IS HUE] seems closely related. 

In focusing on the minimum energy, we obtain something like a phase diagram. Each of the 
constructions used for the upper bound has a different energy scaling law, which can be expressed 
in terms of two nondimensional parameters (see (13.51) 1. By considering which construction is 
best, we divide the parameter space into regions according to the character of the optimal 
construction. Since we obtain scaling laws not prefactors, our analysis leaves some uncertainty 
about the exact locations of the “phase boundaries.” Our lower bound assures that our list of 
constructions is complete, i.e. we have not forgotten any phases. However our results address 
only the energy scaling law - they do not exclude the possibility that some other construction, 
qualitatively different from the ones considered here, could also achieve the optimal scaling in 
some region of parameter space. 

Mathematically speaking, there is nothing novel about using the best of several construc¬ 
tions to obtain an upper bound on the minimum energy. The challenge of finding a lower bound 
that matches (with respect to scaling) the upper bound is however far from routine. Something 
similar has been achieved in model problems motivated by twinning in martensite [23j . uni¬ 
axial ferromagnets mm, the intermediate state of a type-I superconductor [TO], structural 
optimization (26, [27], tension-driven wrinkling in thin elastic sheets (5|, and compression-driven 
wrinkling in thin elastic sheets m ei 0 cei [ 25 ]. All these papers consider nonconvex variational 
problems regularized by higher order singular perturbations, which develop microstructure as 
the coefficient of the regularizing term (call it h) tends to zero. There is a general framework, 
known as “relaxation,” for finding the limiting energy as h —>• o d]. While there is as yet 
no general framework for analyzing the correction to the energy associated with positive h , 
some principles are beginning to emerge (22] ; one is the importance of interpolation inequali¬ 
ties, which play a major role in many of the articles just cited and also in the present work. 
In some settings the relaxed problem is very degenerate, providing little guidance about the 
character of the microstructure. In the setting of this paper (as in 0 ) the relaxed problem has 
a unique solution, which determines the direction (but not the length scale) of the wrinkles. The 
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nondegeneracy of the relaxed problem makes the tension-driven wrinkling problems considered 
here and in [5] quite different from the compression-driven wrinkling problems considered in 

Our work is closely related to contemporary work by the second author and Hoai-Minh 
Nguyen [23] concerning the wrinkling seen in a confined floating sheet [IS]. That problem has 
many similarities to the one considered here: the sheet is in tension, the wrinkling is induced by 
lateral confinement, and the length scale of the wrinkling is smallest at the unconfined edges of 
the sheet. There is an underlying energy, considered already in [18] . The analysis in [23] adopts 
a variational viewpoint similar to that of the present paper, proving matching (with respect 
to scaling) upper and lower bounds on the minimum energy. The construction leading to the 
upper bound uses self-similar coarsening, and is similar to the “type I deformation” discussed 
in Section T6. II The proof of the lower bound shows that changing the length scale costs energy, 
using an argument closely related to this paper’s Lemma [7.81 There is, however, a key difference 
between the hanging drape considered here and the confined floating sheet considered in [23], 
namely: the hanging drape can relieve its confinement by lateral spreading, whereas in the 
floating sheet this is prohibited by the boundary condition. This extra freedom means that 
we must consider spreading as well as coarsening in connection with the upper bound, and it 
means that many new arguments are needed for the lower bound. (The analysis of the floating 
sheet also has complications not present here; in particular, the scale of the wrinkling at the 
unconfined boundary is determined by energy minimization rather than being imposed as a 
boundary condition.) 

It is tempting to think that the energy-minimizing deformations should resemble the ones 
used to prove the upper bound. This is the principle, for example, behind the argument in [3T] 
concerning how the length scale of the folding in a hanging drape varies with height. As noted 
above, however, our rigorous results concern only the energy, not the spatial structure of the 
energy-minimizing pattern. There are in fact relatively few problems involving microstructure 
where the optimal patterns are understood. One is [12] . which considers a model for the 
refinement of martensite twins near an interface with austenite; another is [3], which considers 
a model for the behavior of a sheet near the boundary of a wrinkled region. 

The preceding paragraphs fall far short of a comprehensive review concerning the mechanics 
of thin sheets and the associated mathematical challenges. With respect to wrinkling or folding, 
we have emphasized work in which the length scale of the microstructure varies with position, 
thereby omitting many recent contributions including Gang m m M- Moreover we have 
completely omitted other aspects of thin elastic sheets, such as the formation of localized defects. 
For a broader review concerning the mechanics of thin sheets we refer to [lj. We also mention 
in passing some more mathematical work concerning the energy scaling laws of specific defects 
in thin sheets [a m Eg Eg. 

Returning to our hanging drape, we now offer some heuristics to motivate the analysis that 
begins in Section [2j Here and in the rest of the paper, we often speak of “wrinkles” rather 
than “folds.” This is because we do not expect (and our model does not predict) sharp folds 
similar to creases in a piece of paper; rather, we expect the out-of-plane profile of the drape to 
be smooth (perhaps approximately sinusoidal). Also, we often speak of a “sheet” rather than a 
“drape,” since we model the drape as a thin elastic sheet. 

Recall that the boundary condition at the top involves small-scale wrinkling, which costs 
bending energy. If the thickness of the sheet is large enough to make the bending resistance 
important, then (as noted earlier) one expects to see less bending far from the top. This 
can be achieved by two rather different mechanisms: coarsening of the folds or spreading of 
the sides. A third option - compression of the sheet - is not anticipated since buckling is 
energetically preferred over compression. Our results support the intuition that there should be 
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no compression, since our upper bound uses compression-free constructions while our matching 
lower bound has no such hypothesis. (They do not, however, show that the minimizer is 
compression-free.) 

Gravitational effects oppose both the coarsening of the folds and the lateral spreading of 
the sheet. Indeed, gravity pulls vertically, favoring a configuration that hangs straight down. 
Coarsening increases the amplitude of the out-of-plane displacement, while spreading involves 
horizontal deformation; the presence of either mechanism works against the effects of gravity. 

The deformation of the hanging drape reflects the competition between these effects. A 
gradual deviation from “hanging straight down” is preferred (by the gravitational effects) over 
an abrupt deviation. Therefore we expect the wrinkling to gradually coarsen and the sheet to 
gradually get wider as the distance from the top increases. Our upper bounds are consistent 
with this expectation; moreover they provide guidance (consistent with that in [JI]) concerning 
the rate at which coarsening occurs, and they identify the parameter regime in which spreading 
affects the scaling law. However while the spreading of a real drape is smooth, the spreading 
in our constructions is only piecewise smooth. This is convenient, since it makes it easier to 
estimate the elastic energy; and it is permissible, since our results concern the energy scaling 
law (not the prefactor, and not the character of the optimal deformation). 

The article is organized as follows. Sections [2][5] are all in some sense introductory: Section [2] 
presents our model and discusses in general terms the properties of the low-energy configura¬ 
tions; Section [3] presents our main result - the matching (with respect to scaling) upper and 
lower bounds - and briefly discusses when (in terms of the physical parameters) lateral spread¬ 
ing affects the energy scaling law; Section [4] rescales the energy to decrease the number of 
independent parameters; and Section [3] discusses the “bulk energy,” which plays the role of a 
relaxed problem (its minimum is the limiting value of the energy when h —>• 0). With those 
preliminaries in place, we turn in Section [U] to the upper bound, which is proved by considering 
several candidate deformations. Finally, Section [7] presents our lower bound, which is mathe¬ 
matically speaking the subtlest aspect of the paper. A sketch of the main ideas underlying the 
lower bound is given in Section 17.11 

Notation. We shall denote by C a generic constant, i.e., a constant whose value may change 
throughout the computation. The symbols <, and > indicate that the estimates hold up 
to a finite universal multiplicative constant C, e.g., a < b stands for a < Cb. The tensor 
product u®v is defined as the 3x3 matrix that is component-wise defined by ( u®v)ij = UiVj. 
When f{x,y) is a function, we often use subscripts to denote partial derivatives; for example 
f, x = df/dx and f, xy = d 2 f /dxdy. Finally, when a and b are real numbers, we write a A b for 
the minimum of a and b. 

2 The model 

In this section we discuss the domain, the energy functional, and the boundary condition at the 
top. We also discuss in general terms the expected behavior. 

We assume the drape has (a small) thickness h > 0 and a rectangular shape of width 2W and 
length L (we choose width 2W in order to have a symmetric domain [— W, IF]). We denote the 
domain by := [— IF, W] x [—L, 0]. The curtain is clamped at the top Ty := {(x, y) E fi : y = 0} 
while it is free to move elsewhere. The wrinkles prescribed at will have a (small) wavelength 
u>o; their shape will be specified in a moment (see (12.21) 1. For simplicity we assume that IF = kwo 
for some (typically large) integer k. 

As usual in elasticity, the stable configurations of the drape are local minima of an “energy 
functional,” obtained by adding the elastic energy and the work done by gravity. Our goal (as 
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discussed in the Introduction) is to understand how the minimum energy scales with respect to 
the physical parameters. Evidently, we are studying the energy of the ground state. 

For the elastic energy we use a geometrically linear Foppl-von Karman model, and we take 
Poisson’s ratio to be zero. This is, admittedly, a qualitatively accurate model not a quantita¬ 
tively accurate one: for real materials Poisson’s ratio is usually not zero, and the Foppl-von 
Karman framework is only appropriate when the out-of-plane deformations have small slope. 
We believe, however, that our choice captures the essential physics of the phenomena we wish to 
study. This view is supported by the mechanics and physics literature on wrinkling, where the 
Foppl-von Karman framework is widely used. It is also supported by the mathematics literature 
on thin elastic sheets, where energy scaling laws initially derived using a Foppl-von Karman 
model have been shown to hold also in more nonlinear settings, see e.g. [2 [13]. (For further 
discussion about the appropriateness of the Foppl-von Karman framework see e.g. [32].) Based 
on the preceding considerations, our energy functional is: 

Eh(u,Q= [f \e{u) + V£| 2 + /i 2 |V 2 £| 2 dxdy + T ff u y dxdy, (2.1) 

JJn z JJn 


where u = (u x , u y ) denotes the in-plane displacement, £ is the out-of-plane displacement, e{u) = 
\/u+\/u T - g |j ie symmetric gradient of u, and r > 0 is a given parameter (the ratio between the 
gravitational constant and Young’s modulus of the elastic material). For our geometrically 
linear Foppl-von Karman model to be reasonable we want the curvature to be much smaller 
than 1/h. This means that any length scale in the deformation (in particular the period wq of 
the prescribed wrinkling at the top) should be larger than h. 

We turn now to the boundary condition u(x, 0) and £(x,0) imposed at Ty (the top of the 
sheet). The out-of-plane deformation £ should be periodic with period wq (note that since 
we assume W = kwo, Ty is filled by exactly k periods). Moreover we want to avoid strain 
in the horizontal direction, and the deformation should achieve a specified overall horizontal 
compression factor A. Finally, it is natural to choose u(x, 0) and £(x,0) so that the bending 
energy is minimized subject to these constraints. Focusing initially on the two-period interval 
[— wq, wo], we seek rto,£o : [—'ixo,ico] —> M such that 

uo,x( x ) +£o,x( a 0/2 = o, 
u 0 (-w 0 ) - u 0 (w 0 ) = 2Aw 0 , 


and such that the bending energy 


h; 


/ WO 

' 

-WQ 


£o,xx( x ) dx 


is minimized subject to these constraints. Using the method of Lagrange multipliers one finds 
that the choice 

£o(x) = “ 0v/ ^ sin( 27 rx/w 0 ) 


7r 

is optimal. Extending £o by periodicity, we are led to impose the boundary condition that 

£(x,0) := sm(2nxwQ 1 ), u x (x, 0 ) := f \^ x (t, 0)| 2 dt, u y (x, 0 ) = 0 (2.2) 

71 1 Jo 

at Tt- 

The elastic energy of our boundary condition is h 2 f U w £ 2 XX dx = 2Wh? A{Stt 2 Wq 2 ). For the 

trivial planar deformation £(x,0) = 0, u x (x,0) = —Ax the elastic energy f W w r u 2 x dx is 2IUA 2 . 
Our boundary condition has lower energy than the trivial one when 


h < u;oV / A(87r 2 ) 1 / 2 . 


(2.3) 
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Thus it is reasonable to prescribe (12.21) only if (12.31) is satisfied. 

2.1 The expected behavior 

We now discuss the expected form of a deformation with small energy. This is, in effect, a 
description of the deformation associated with our upper bound - which achieves the optimal 
scaling, according to our lower bound. 

Since we do not prescribe boundary conditions on the lateral part of the boundary, the sheet 
is free to get wider (and, therefore, to relax part of the confinement forced by the boundary 
conditions at the top). Using two different constructions, we will show that the energy required 
to significantly change the value of u x (i.e. to release the lateral confinement) over the length Iq 
from the top of the sheet is at most of order W 2 A 2 min(Wi^" 1 , IT 3 /^ 3 )- The energy required to 
significantly decrease the amplitude of the out-of-plane displacement £ scales like W(Awq)tLIq . 
If the sheet releases the confinement, then both of these terms contribute to the energy. 

Assuming the sheet avoids compression, besides getting wider towards the bottom it needs 
to waste some arc length in the horizontal direction in the region where it is confined. This 
can be done by buckling out of the plane in wrinkles with an “average period” w = w(y). The 
boundary condition at Tt sets w(0) = wq. The bending energy of a deformation consisting of 
wrinkles with period w is 16ir 2 WAh 2 w~ 2 . It is obvious that this term prefers to increase the 
period w as fast as possible towards the bottom of the sheet (away from T^), but a rapid change 
in the length scale w would require a large change in the amplitude of the wrinkles (large £ y ). 
Since the sheet is stretched in the vertical direction by gravity, the part of the energy coming 
from stretching in the vertical direction prefers small £.,y. In the end, the competition between 
these two preferences determines the rate at which the length scale w(y) increases. 

In our constructions, the variation of w(y) is achieved using “building blocks.” Each building 
block is a deformation defined on a rectangle of width w with sinusoidal profiles at the top and 
bottom boundaries with period w and 3 w respectively. The idea of using such building blocks 
was already present in [SUSHIS]; our building blocks are called wrinklons in m- As we’ll see 
in Section [S energy minimization requires the height l and width w of a building block to be 
related by l ~ w 2 \frLhr x . 

The optimal number of building blocks - i.e. the number of generations of coarsening - 
depends on several parameters. It is an increasing function of L (in a longer sheet there is room 
for more generations of coarsening), and an increasing function of h (a thicker sheet is harder to 
bend, so it decreases the amount of bending faster). On the other hand, the number of building 
blocks is a decreasing function of wq (finer wrinkles require more bending, so a sheet with finer 
wrinkles at the top prefers to coarsen the wrinkles faster), and a decreasing function of r (in a 
heavier drape the effect of gravity is stronger, so out-of-plane displacement is more expensive, 
which implies that coarsening is also more expensive). 

If the sheet is very long, the coarsening process may finish at a height above the bottom of 
the drape. As one goes further down, the out-of-plane displacement becomes an affine function 
with the correct slope. From this point on, there are no more contributions to the energy and 
the sheet does not change its shape (except for the vertical deformation due to stretching). 

Our discussion has emphasized the properties of the deformation in the horizontal direction. 
In the vertical direction the situation is much simpler. Since the sheet is pulled down by gravity 
independently of the horizontal position, we expect the vertical deformation to be independent 
of x (up to a small correction due to the wrinkling). In Section [5] we will formulate a one¬ 
dimensional variational problem which will be used to identify the main part of the optimal 
vertical displacement u y . 
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3 The main result 


This section presents our main result, and provides some discussion to help elucidate its conse¬ 
quences. 

Theorem 1. Assume 


tL > 4, (3.1) 

hA~ l/2 < w 0 < (2c w )W, (3.2) 

where c w > 0 is a small universal constant. Then there exist universal constants Cub > Clb > 0 
such that for any deformation («,£) which satisfies (12.211 we have 

- 1- T 2 L 3 (2W) + C LB e < min E h (u,Q < ~ t 2 L 3 (2W ) + C UB e (3.3) 

1Z (u,£) J .Z 

where 


e := W A min 


< hy/rL log ( 


hyfrL log ( 


( — 2 ( hL 




Wn 




A 4W 2 + 1 , 


mm 

l£(0,L) I 


hi 


\wqVtL 


+ 1 ] + WqtLI 1 + WA min ( ( -|-j 



(3.4) 


Remark. The first inequality in m is motivated by (12.311 (dropping the constant (8-k 2 ) x ! 2 
compared to (12.31) does not change the scaling law). The second inequality in (13.21) (with c w <C 1) 
says that the period of the wrinkling prescribed at the top is much smaller than the width of the 
sheet. 

We see that the minimum of the energy in Theorem [7] consists of two parts. One is the 
‘ l bulk energy” —-^t 2 L 3 (2W), which comes (as we’ll show in Section [5J) from the stretching of 
the sheet in the vertical direction on account of gravity. The other is the “excess energy” due 
to positive h; it is of order e. We observe that e does indeed vanish as h —>• 0/ this is consistent 
with the fact that wrinkling uses less energy in a sheet of smaller thickness. 


Remark. The paper J22F distinguishes between “heavy sheets” and “light sheets,” and shows 
that the rate at which wrinkles coarsen is different in the two cases. The assumption means 
that our drapes are “heavy sheets. ” 


The excess energy e is the minimum of two different terms - the first is the energy of the 
construction when the sheet does not (significantly) release the lateral confinement; the latter 
is the energy of the construction when the sheet releases (much of) its lateral confinement. It 
is natural to try to understand, in terms of the physical parameters, the regimes in which one 
or the other term dominates. To this end we introduce the non-dimensional parameters: 


hL 

Wq\/tL' 


P := 


r := 


l_ 

T 


We focus on the case when the sheet is not extremely long, in the sense that a = < 

Then we get the following formula for the “average” excess energy: 


e 


LWA 

= min aB 2 log (a + 1), min 
\ re(0,l) 


a/3 2 (log(ar + 1) + (ar) A min 
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We first neglect the last term A min ^ (j^) , (jy:) 3 ) and investigate when 

min logfar + 1) + (or) -1 < log (a + 1). (3.6) 

re(0,l) 

Since r € (0,1), for (13.61) to hold, the function J-(t) = log(f +1) + t~ 1 must achieve its minimum 
at a value f m j n that is smaller than a. Therefore (|3.6I) can be satisfied only if 

l + \/5 

Since the last term in (13.51) (the one we neglected) is positive, it is clear that when a < (l+\/5)/2 
the first term in (13.41) is larger, i.e. the sheet does not prefer to release most of the lateral 
confinement. On the other hand, if a is considerably larger than (l+\/5)/2 we expect the second 
part of the right hand side in (13.51) to be smaller provided the neglected term is sufficiently small, 
in which case most of the lateral confinement should be released. 

4 The rescaling 

The scaling law (13.31) depends on numerous parameters: the sheet’s length L, its width W, 
and its thickness h , the “gravitational” coefficient r, and the period of the wrinkling at the 
top wq. To simplify the analysis, it is convenient to reduce the number of parameters by 
nondimensionalizing the problem and by using its special structure to eliminate A. 

The effect of nondimensionalization is that it permits us to consider only the case W = 1/2. 
This is achieved by measuring length in units of width. Explaining in detail: if ( u ,£) is a 
deformation defined in (—IT, IT) x (—L,0), we define a new deformation (v,y) by: 

v(x, y) := {2W)~ 1 u(2Wx, 2Wy), y{x, y) := {2W)~ 1 ^2Wx, 2Wy). 

The deformation (u,/i) is defined in (—1/2,1/2) x (— L/(2W), 0), and we have 

-^'h,l/2,L,f,wo ( V ’ = (^W) Eh,W,L,T,u> o{ u iC)i (4-1) 

where we have listed all relevant parameters for the energy as indices. The rescaled parameters 
are: 

h = (2W)~ 1 h, L = (2W)~ 1 L, t = (2W)t,w 0 = (2W)" 1 w 0 . 

This nondimensionalization replaces L, h, and wq (which are lengths) by their quotients with 
2 W ; similarly, it replaces r (which has dimension length -1 ) by its product with 2 W. We observe 
that (14.11) is consistent with (13.31) . (13.41) . Therefore to prove Theorem [1] it is sufficient to consider 
the case W = 1/2. 

Our second reduction uses the special structure of our energy functional - specifically, the 
fact that the membrane term is quadratic in u and quartic in £, while the bending term is 
quadratic in £, and the gravitational term is linear in u — to eliminate the parameter A. Defining 

v(x, y) ■■= A -1 ri(x, y), n(x, y) := A -1/2 ^(x, y), 

we see that 

^h,T ,1 (^’ A 4 ) = ^ ^h,T, a{ u iC)i 

where h = A= A -1 t. This relation is again consistent with Theorem [1] Therefore in 
proving the theorem it is sufficient to consider A = 1. 

In the rest of the paper we will assume that A = 1 and W = 1/2. Note that in this case the 
second hypothesis (|3.2|) of Theorem [T] becomes 

h < vjq < c w . (4-2) 
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5 The bulk energy 


The first term in our energy scaling law (13.3D is the limiting value of the minimum energy as 
/i —>- 0. While the proof of this assertion is given later on (in the course of establishing our 
upper and lower bounds), the present section lays necessary groundwork by considering a “bulk 
energy” that includes only vertical stretching and gravity. 

To motivate the definition of the bulk energy, we begin by substituting h = 0 in the definition 
of our functional Eh- Gravity pulls downward, so we expect the vertical displacement to satisfy 
Uy t y > 0, and assuming this we have 


Eo(u,0> // \u y , y (x,y)+(, 2 y (x,y)/2\ 2 dxdy + r u y (x,y)dxdy 
JJu JJ n 


> jI \uy >y \ (x,y) +Tu y (x,y)dxdy. 


Therefore it is natural to consider the bulk energy functional, defined by 




+ rfdx dy, 


(5.1) 


where / : [—L, 0] —> R, /(0) = 0 (/ plays the role of u y ). It is easy to find the unique minimizer 
of B: 

f{y) = (ry 2 + 2rLy)/ 4, min B = B(f) = t 2 L 3 . (5.2) 
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The following estimate for f <y (y) = r(y + L )/2 will be useful later: 


f,y(v ) ^ rL / 4 i f V € [ -L/2 ,0], 

0 < f, y (y) < rL/2 iiye [-L,0]. 

We see that to have f >y (y) > 1 for y > —L/2 we assumed ()3.1I) . For u y (x,y) := f {y),u x (x,y) := 
u x (x,0), £(x,y) := £(x,0) we have 


Eo(u,t) — I T Uy )X + ^x^,y\ /2 + l^y,2/ T ^ y /2\ + Tu y dx dy 


\ u y, y \ 2 + ru y dx dy = 


thus min / u ^Eq(u,^) (subject to u Vty > 0) is attained and is equal to B(f). We will see later 
that for a deformation (it, £) to nearly minimize Eh, the vertical part of the deformation u y 
must be close to /. 


6 The upper bound 

This section proves our upper bound for the minimum energy. As shown in Section [U it suffices 
to consider W = 1 and A = 1. Our goal is therefore to show that 

min(it, £) < + C UB e, (6.1) 

provided (|3.1D and (14.2D hold, (it,£) satisfy the prescribed boundary condition (12.2D . and e is 
defined by (13. 4D . 

The proof of (16.ID uses three types of deformations. Our type I deformation involves self¬ 
similar coarsening of the wrinkles, while the sheet remains horizontally confined. The coarsening 
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is achieved through multiple generations of building blocks; the period of the wrinkles is tripled 
in each generation. If the vertical length of each generation is chosen in the optimal way, its 
contribution to the energy is of order hy/rL. If the sheet is very long - more specifically, if 
we reach a generation of wrinkles with period comparable to the width of the sheet - then 
we change (within one generation) the out-of-plane displacement to an affine function (since 
when £ is an affine function with the correct slope all the terms contributing to the excess 
energy vanish). This construction involves no lateral spreading: the values of u x (— 1/2, y) and 
u x (l/2,y) are independent of y. The vertical deformation u y agrees with the minimizer / of 
the bulk energy B. 

To be energetically efficient, the coarsening process associated with a type I deformation 
needs some room. We discuss the case when L is too small for coarsening at the end of Section 
IMl Since there is no room for coarsening, in this setting the type I deformation keeps the 
profile of the wrinkling independent of y. 

Our second type of deformation is a modification of the first. It also involves self-similar 
coarsening of the wrinkles, but in contrast with the first type the horizontal confinement is 
relaxed at some point. To describe it, consider a type I deformation, and choose a particular 
generation in the coarsening process. We denote by n the order of this generation (i.e. there are 
n — 1 generations above) and by l n the vertical length of the building blocks in this generation. 
Our type II deformation is identical with the type I deformation through the first n — 1 gen¬ 
erations of the coarsening process, but different starting at the n-th generation; it completely 
releases the horizontal confinement of the sheet within the n-th generation (the value of u x at 
the extremes ±1/2 changes from ±1/2 at the top of the nth generation to 0 at the bottom). 

Full details of the type II deformation are given in Section RT21 but here is a sketch. If we keep 
u y = /, changing u x from order 1 to 0 over the length l n results in the term \u XiV ± u V)X ± £, x £,y| 2 
being of order l~ 2 (since u VtX = 0 and the term involving derivatives of £ is not larger than 
u Xj y). Integrating this term over the domain of size l n , we find that this mechanism for releasing 
the horizontal confinement has an energetic cost of order 0(l n * (l/l n ) 2 ) = O(In 1 )- A different 
possibility is to set u VtX := —u x>y — £ >x £ >y . This obviously eliminates the term \u x , y + u ytX + £ tX €,y\ 
from the energy, but it increases the term \uy iV ± £^/2|. For this alternative mechanism of 
releasing the horizontal confinement the energetic cost turns out to be 0(l~ 3 ). 

Our type III deformation is similar to a type II deformation with n = 0. It flattens the 
sheet and eliminates the horizontal confinement over an interval of height starting near y = 0. 
As with a type II deformation, there are two versions of this construction, corresponding to 
different choices about which membrane term should be made to vanish. 

6.1 The construction of a building block and the type I deformation 

In this section we construct a deformation consisting of several generations of building blocks, 
with the length scale of wrinkling being tripled in every generation. We first define a building 
block, which is a deformation (v, fi) (the horizontal and the out-of-plane displacement) defined 
on [0,1] x [0,1] with wrinkles with period 1/3 at the top and with period 1 at the bottom with 
periodic lateral boundary conditions (see | 30] for a numerical study of the optimal shape for one 
such building block). Our Type I deformation is then obtained by patching together rescaled 
versions of the building block (the boundary conditions for the building block were chosen to 
make this possible). 

The building block is a deformation (v, n) : [0,1] x [0,1] —>• M 2 . It has a sinusoidal profile of 
period 1/3 in a neighborhood of y = 0, and a sinusoidal profile of period 1 in a neighborhood 
of y = 1. There is neither compression nor tension in the x-direction, and both the membrane 
and bending energy are finite. Saying the same in mathematical terms: both x ± v(x, y) and 
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fi(x, y) should be 1-periodic in x, and they must satisfy: 


Mu y ) : 

= -^sin(67rx), w(x,y) = -| 

[ dt, 

0 

[0,1], 5 € [0,1/4], 

(6.2) 

Muy) 

1 1 1 

= -sin(2vrx), u(z,y) = -- 1 

7r 2 J, 

Im,®(^, i)I 2 dt, 

0 

[0,1 ],ye [3/4,1], 

(6.3) 


v, x (x,y) + -\fi :X (x,y)\ 2 = 0 for (x,y) € [0,1] 2 , 

M < 1, H < 1 

(6.4) 

E'm ■ — 

// |u y + //,x//,y| 2 + \n,y\ 2 dx dy < oo, E b := 

7( o,i) 2 

// V 2 /i 2 dxdy < 00 . 

//(0,1) 2 

(6.5) 


Such a building block can be obtained as follows: let 51,52 : [0,1] —>• [0,1] be smooth functions 
which satisfy 

5 i(y) = l, y €[0,1/4], 

52 ( 5 ) = 1, y €[3/4,1], 

9i(y) + gl{y) = i, ye[o, 1 ], (6.7) 

Wi(y)\ +3|52(y)l < 3tt, y € [0,1]. (6.8) 

We define 

fj>(x,y) ■= ^^sin(67rx) + ^^sin(27rx), v(x,y) := [ \y, x (t, y)\ 2 dt, (x, y) <E [0, l] 2 . 

ovr 7r 2 J o 

Using the properties of 51 and 52 it is easy to verify (16.2116.511 . Indeed, (16.211 and (16.311 follow 
from (16.611 and (16.41) follows from (16.71) and (16.81) . Finally, (16.51) is a direct consequence of the 
smoothness of 51 and 52 . 

To define the type I deformation we patch together rescaled versions of the building block, 
and use u y (x,y) := f(y) for the horizontal displacement. The new rescaled displacement u x ,£ 
(defined on a rectangle of size w x l with the upper left corner at (xq,5o)) is defined as: 


u x (x,y) := x 0 + wv 


x-x 0 y - 50 


w 


l 


Z(x,y) ■= wy 


x-xo y-y 0 


w 


l 


From \y\ < 1 we get |£| < w. The period of the wrinkles at the bottom and upper end of this 
building block are w and w/3, respectively. We also have 


u. 


:,x + £,®/2 — v, x + lf x /2 — 0, 


w 


W' x,y T VT y )X T ^x^,y — ^ (V,y ~f" Ah^/hl/) i 
{u y ,y + i 2 y / 2 ) 2 - \f,y \ 2 = f,y$y + ^/ 4 = {w/lf (f, y ff y + {w / if ^ / 4) , 


f.xx — VC y,xxt ^,xy — l l-t'jXyi ^,yy — t^,yy 

In what follows we will always keep w < l. Then (16.41) implies \{w/if yf y /4\ < ff /4 and 
|V 2 e| 2 < uU 2 |V 2 /v| 2 , and the elastic energy of one building block (over the rectangle of size 
w x l) is bounded by 


-1 . 


;- 2 . 


2 wl 


(j) 2 (i + /')^m + h 2 W~ 2 E b 


+ II (f ' 2 + rf). 


(6.9) 


When summed over the whole domain U the second term becomes exactly m- We are 
interested in the first term which is the excess energy due to positive h (roughly speaking: the 
energy due to wrinkling). 
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We will patch together several generations of building blocks. Let N be the number of 
generations (to be chosen later), and let l n and w n := 3 n wo for n = 1,... ,N be the length 
and width of the building block in the n-th generation, respectively. We may assume without 
loss of generality that 3 K wo = 1/2 for some integer /Ljj Since one block has width w n and the 
width of the sheet is 1, each generation has w~ * * 3 identical building blocks. We sum the first 
term in (EHD over all building blocks to get 




(1 + f')E m + h 2 w~ 2 E b 


N 


n =1 


—r-TLE m + h 2 w n ~l n E b , 


( 6 . 10 ) 


where we used (15.31) and (13.11) to obtain 1 + /' < 1 + tL/2 < tL. 

We know that for (16.101) to be small we want the two terms on the right-hand side to be 
2 

of similar value, i.e. j 21 TLE m ~ h 2 w~ 2 l n E b . From w n = 3 n u;o (the period is tripled in each 


generation) we obtain l n ~ 


9 n w 2 0 



Motivated by this we set 


l n := w^VrLh 1 = ^w^y/rLh 1 . 


( 6 . 11 ) 


Using m we can now verify the previously used assumption w n < l n : 

l n /w n = 3 n VrLwo/h > VtL > 2. (6.12) 


Let us first assume that the sheet is not very long in the sense that 

I< 


in > l. 


(6.13) 


n= 1 


Then the length of the sheet L can be expressed as the sum of the lengths of all generations 
of building blocks, i.e. Yln =l ^ ~ L. Let us dehne N to be the smallest integer such that 
J2n =i > L. It follows from (16.131) that N < K and so w n < 1 for n = 1,..., IV. Using (16.111) 
we obtain 


N = 

where we assumed that 


!ogg 


8 hL 


9 w%VtL 


1 


log 


hL 


\VtL 


+ 1 


w; 


(6.14) 


h > WqVtL/L. (6.15) 

From the definition of l n we compute that one generation of wrinkles (re” 1 identical build¬ 
ing blocks) costs Chy/rL. Therefore, for L not too large the energy (16.101) is bounded by 
+1^ and 


ChyfrLlog ' hL 


VtL 


min E h (u,£) < --}-t 2 L 3 + Chs/rL log [ + 1 ) . (6.16) 

!2 yw^y/rL J 

Now we treat the case when (16.131) is false, i.e., Lq := Yln=i ^ < L. For y € [—Lo>0] we 
define the deformation the same way as before. For that we need K generations, where each 
contributes to the excess energy by a multiple of hyfrL. Since K is defined through 3 K wo = 1/2, 

Tndeed, if (before nondimensionalization) there is a I\ such that 3 K ^ 1 wo < W/2 < 3 K wo, then we can 

consider a slightly wider drape, of width W such that W /2 = S K wo (for an upper bound it does no harm to 
increase the domain). After nondimensionalization (which divides all lengths by IF), we get a problem in which 

3 k w 0 = 1/2. 
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the excess energy for this part will be of order h^/rL log^g 1 ). To define the deformation for 
y G [-L, —Lp], we first observe that at y = —Lp we have 


£0, -L 0 ) = —sin(2vrx), u x (x,-L 0 ) = ~- £ 2 x (t,-L 0 ) dt, zG[0,l]. 


7 r 


(6.17) 


We set Ik+i '■= VrL/h, and for x G [0,1 ],y G [—L,Lq) we define 


£{x,y) = 9i 


y + L 0 \ 
l K +1 J 


sin(27rx) 
7r 



y + L 0 \ 

h<+\ ) 


V2x, 


u x (x,-L 0 ) = 


where 51,32 are functions 31,32 from (16.6ti6.8D . extended respectively by 
With the above choice of Ik +1 it is easy to compute that the contribution 
in [—(To + Ik+i), —Lp] is of order h\frL , while for y < —(To + Ik+ 1 ) we 
u x {x,y) = x, so there is absolutely no contribution to the excess energy, 
larger than J2n=i ^n, the contribution to the excess energy is at most 


\ J Z%(t,~Lo)dt, 

(6-18) 

0 and 1 into [l,oo). 
to the excess energy 
have £(x, y) = \/2x, 
We see that if L is 



Hence, together with (16.161) we get that 


min E h (u,£) < 


- ^r 2 L 3 + ChVrL log ^u; 0 2 = A l^j 



(6.19) 


( 6 . 20 ) 


If (16.151) does not hold, i.e. if h < Wq\VtL/L, it should be better energetically to just 
propagate the deformation prescribed at Tip. In this case we set 


u x (x,y) := u X! q(x) = u x (x,0), u y (x,y) := f(y), £(x, y) := £ 0 (x) = £(x, 0), (6.21) 


and we get the total energy bounded by 


min E h (u,£) < 


1 


- — T L + Ch w o L 


( 6 . 22 ) 


Finally, since log(l + t) > t /4 for t G (0,1) and < 1 follows from h < WqVtL/L, we see 
that (16.201) and (16.221) can be rephrased as 


min E h (u , £) < - ^t 2 L 3 + ChVrL log ^rc 0 2 A 1 j + 1 j ■ (6.23) 


Remark 6.1. Our type I deformations are equivalent to the ones discussed for heavy sheets 
in m ■ In particular, our relation (16.111) is the equivalent in our notation of equation (4) in 
m, giving the optimal “length of the wrinklon.” Since {l n } is a geometric series (with ratio 
greater than one) we have Ip + • • • + l n ~ l n , so length scale of the wrinkles at height y - 
call it w(y ) - can be read off from (16.111) : it satisfies \y\ ~ w 2 (y)VrLh 1 . Rewriting this as 
- ~ (rL) -1 / 4 ^) 1 / 2 we see that it is the equivalent in our notation of equation (5) in fJTf. 


6.2 The type II deformation 

This construction is a modification of the previous one; besides coarsening the wrinkles it also 
releases the horizontal confinement of the sheet. Let ( u , £) be the deformation from the previous 
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section, constructed using N generations of coarsening for some N < K. We choose a particular 
n such that 1 < n < N. We denote the y-coordinate of the k -th building block by 

fc-i 

$k ■— ^ j li- 
i= 1 

For y E [— s n ,0] the new deformation (17, C) coincides with («,£): 

U(x,y) := u(x,y), 

C(x,y) ■= £(x,y). 


For y E [—L, — s n ] we define 


U x (x,y) := ip 2 - 


s n + y 

In. 


u x (x,-s n ), C (x,y):=<p[- 


s n + y 
In 


£(x,-s n ), 


where ip : [0, oo) —>• [0,1] is a smooth decreasing function which satisfies \<p'\ < 2 and 

y(t) = 


1 t E (0,1/3), 
0 t E (1, oo). 


Then for y < — s n we have 

U x , x (x,y) + ( 2 x (x,y)/2 = p> 2 ^_£ZLI tl'j ■ (u x>x ( X ,-s n ) + gf x ( X ,-s n )/ 2) =0. 

So far we did not define U y for E (— L, —s n ). We will do this in two different ways, thereby 
obtaining two different upper bounds for the energy. The first way simply sets 


U y (x, y) := u y (x, y) = f(y). (6.24) 

For y E [—s n ,0] the excess energy is estimated by the second term in (|6.23l) (with L replaced 
by s n in the numerator), and so we just need to estimate the excess energy in the part of the 
domain y < —s n . We see 


U x ,y T Uy tX + C,xC, 2 /l — | (2ip<p )(( S n y)/l n )u x (x, S n )l n 

+ (</V)((- s n - y)/ln)ln 1 ix(x, S n )£(x, -S n )| 2 < Cl~ 2 , (6.25) 


where we used that |£(-, — s n )|, |£ )X (-, — s n )|, \u x (-, — s n )|, |</9 / |,|</2| < C. Using (16.121) we see that 
|C,2/1 < l^'IICIAn < 2 w n /l n < 1. So ([OP implies 


I Uy,y + C,y/2| 2 - \f,y\ 2 = /, 2 /C 2 + Ci/4 < rL\C y \ 2 = TL\<p'\% 2 e(x, s n ) < 4 tL1^ 2 (x, -a n ). 


Since |£(x, — s n )| < w n and l n = w 2 \frLh 1 , the previous inequality implies 

2 

\Uy,y + C 2 /2| 2 - \f,y\ 2 < 4tL/- 2 C 2 (x, -s n ) < 4tL^I~ 1 = 4 hWZ/l n . (6.26) 

It remains to estimate the bending energy. For y E \—L, s n ] we have by (16.121) 

\C,xx( x , y)\ = ip \^xx(x, s n)/ 1$ w n = vVT/( hl n ), 

\(,xy(x,y)\ 2 = tp ,2 \ix(x,-s n )\ 2 l-' 2 < l~ 2 < w~ 2 = VrL/(hl n ), (6.27) 

1C, yy(x,y)\ 2 = <p" 2 Ux,-s n )\% A < w 2 J~ A < w~ 2 = VrL/(hl n ), 
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and so 


h: 


r f 

J maxf—L,—Sn-n ) •' — 


' max(—L,—s n _|_i) J-1/2 

Now we combine (16.281) with (16.251) and (16.261) to obtain 


1/2 \/tL 

V 2 C| 2 dxd y < h 2 (s n - s n+ i)^^- = hVrL. 


hl r 


(6.28) 


E h (U,()<--^r 2 L 3 + C (hV^L 




log 


hs r 


W, 


\JtL 


+ 1+1 


In * l n 


-2 


< --^t 2 L 3 + C ( hVrLl 




og 


Hlr 


+ 11+C 1 ), (6-29) 


WtL 


w. 


where we used that s n and l n are comparable to replace s n with l n in the logarithm, and the 
fact that hl n l(w\\frV) = 9 n > 1. 

The alternative way to choose U y will give the different upper bound 


E h (U,{)<-—T 2 L 3 + C 


The essential idea is to choose U y so that 


( hV tL L 


hlr 


og 


vyk\[rL 


+ 1 I + Z, 


-3 


Ex,y + E y ^ x + C ,xC,y — 6 


(6.30) 


(6.31) 


for y < —s n . Whereas our previous choice incurred a substantial energetic cost from the term 
\U x ,y + Uy !X +^ x C , y \ 2 , our alternative choice makes this term vanish, at the expense of an increase 

in \ U y,y + C,y/2| 2 - 

We can assume that l n > 1, since otherwise (16.301) is worse than (16.291) . To satisfy (16.311) we 
define 

rx 

U y (x,y):=f(y)~ U x , y (s, y) + ( iX (s, y)(, y (s, y) ds (6.32) 

Jo 

for x € (—1/2,1/2) and y < —s n . Then 

2 f-Sn + y 


Uy,y{x,y) = f, y {y) - V 


\Uy,y{x,V)\ < \f,y(y)\ + 


lr 


2 f -Sn + y 
In 


,yy 


/ :X 

u x (s, y) + £, x (s, y)£{s, y)/2 ds 

rx 

/ u x (s,y) + ^ x (s,y)£(s,y)/2ds 
Jo 


,yy 

Y/7-2 


< tL/2 + C\(ip ), yy \ < tL/ 2 + Cl n < tL, 
where we used that l n > 1 and (15.31) . Therefore 

\U y , y (x,y) + CK x ,y)/2\ 2 - \f, y \ 2 = \U y , y (x,y)\ 2 - \f, y \ 2 + U y , y {x,y)Q 2 y (x,y) + ^ y (x,y)/A 

< \Uy,y(x,y)\ 2 - \f, y \ 2 + CtL( 2 v . 

The contribution to the energy from CrL(^ 2 y can be estimated as in (|6.26p : it is at most Ch\frL. 
Since / is the minimizer of the bulk energy B we have that 


It/, 


y,y\ 


+ tU v dx dy ) — 


\f, y \ 2 + rfdxdy 


= B(Uy) - B(f) = - ( D 2 B(f)(Uy - /), Uy-f)= || Uy,y ~ f,y \ | | 2 („) 

Sn + y 


< 




,yy 


dxdy < l n * l n A = l n s . (6.33) 
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The estimate for the bending energy (I6.28|) remains valid, and so the combination of (16.331) 
with (16.311) implies (16.301) . 


6.3 The type III deformation 

The construction in this section is closely related to the discussion at the end of Section 16.11 
which considered the consequences of “propagating the deformation prescribed at Ft” (see 
(16.21 1) 1. Here we do something similar, but we release the horizontal confinement by a mechanism 
similar to that of Section 16.21 

Choose any l E (wq,L), and consider 

U x (x, y ) := y 2 (~y/l) u x (x, 0), U y (x, y) := f(y), ((x, y) := ip (- y/l ) £(x, 0). 

Then \U x>x (x,y) + (, 2 x (x,y)/ 2| = 0 for all (x,y) € H, and 



\U x ,y + Uy !X + C,xC,y\ Z dxdy<l*l 2 = l X , 


r 1/2 

\U y ,y + C,y/2| 2 — \f,y\ 2 dx dy < tLI 1 |^(x,0)| 2 d x<wlrLl X . 

J- 1/2 


Since l > wq, a calculation similar to (16.271) shows that h 2 ff |V 2 C| 2 < h 2 w 0 2 l. Combining these 
estimates gives 

E h (U, C) < - Y -T 2 L 3 + C [h 2 wfl + r 1 + wlrLl - 1 ) . (6.34) 

In Section [6.21 we considered two different ways of extending U y . The preceding calculation 
is like the first, but we can also consider the second. Using (16.321) to define U y and proceeding 
as above we find the estimate 

E h (U, C) < - j^r 2 L 3 + C (h 2 WQ 2 l + r 3 + wlrLl- 1 ) . (6.35) 

Since the right-hand side of (16.341) as a function of l E (0, L) attains its minimum for l > wq, 
we immediately observe that (16.341) holds for all l E (0, L). Finally, for the same reason (16.351) 
holds for all l E (0, L ) as well. 

Taken together, the upper bounds (16.231) . (16.291) . (16.301) . (16.341) . and (16.351) establish (16.11) . 
Thus we have proved the upper bound half of Theorem |T] 


7 The lower bound 


This section proves our lower bound for the minimum energy. As shown in Section [4] it suffices 
to consider A = 1 and W = 1. Our task is therefore to show that if m and (14.21) hold and 
(ii,£) satisfies the prescribed boundary condition (12.21) . then the excess energy satisfies a lower 
bound of the form 


6 := E h (u, 0 - B(f) > C LB min (hVrL log (w 0 2 = A 1 j + 1 j , 


mm 


i€(0,L) I 


< hVrL log ( 


hi 


\wqVtL 


+ 1 ] + WqtLI 1 + min (l l , l 3 ) 


Here / is the minimizer of the bulk energy functional B , and Clb > 0 is a (sufficiently small) 
constant that does not depend on the parameters of our problem (several smallness conditions 
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on Clb will emerge in the course of the proof). We will argue by contradiction; in fact, our 
strategy is to assume that 




<5 < Clb min hVrL log 


( hL 




^Wq^tL 


+ 1 


mm 
le(o,L) I 


< hVrL log 


( hi 


\^WqVtL 


+ 1 J +WqtLI 1 + min (l 1 ,1 3 ) 


(7.1) 


and to prove using this smallness condition on <5 that our lower bound (the opposite inequality) 
must hold. 


Remark 7.1. We can assume without loss of generality that the deformation (tt,£) is smooth, 
since mollification has only a small effect on the energy. (It is important here that our goal is 
the scaling law, not the optimal value of the prefactor Clb-) 


7.1 The idea of the proof 

Before beginning the proof in full detail let us sketch the main steps. 

• Since the sheet is stretched vertically, it prefers not to change its out-of-plane displacement 
£. The situation is similar to a stretched rubber band, whose preferred configuration is the 
straight line joining its endpoints (and for which deviation from this configuration costs 
additional elastic energy). We also know that when h = 0 it is optimal to have u y = /, 
and we expect similar behavior for h > 0 (modulo small adjustments due to wrinkling). 
We make these ideas quantitative in Lemma 17.21 

• In some cases we expect the sheet to spread laterally, releasing the horizontal confinement 
prescribed at Ty. Mathematically speaking, spreading entails decreasing the value of \u x \, 
which is of order 1 at y = 0. If the value of u x is decreased significantly over a length Iq, we 
expect u XjV to be of order Iff 1 . The term u XtV appears in the energy in \u x ^ + Uy tX + f x f y \. 
So we are left with two alternatives: either u yjX and £, x £,y are negligible, and u X}V being 
of order If 1 over a domain of size Iq makes the excess energy at least of order If 1 ; or else 
one of the terms u V:X or f x ^ y is of order If 1 . In the second case other terms in the energy 
come into play. In Lemma 17.31 we show that if the sheet releases a non-trivial part of its 
confinement over length Iq, then the excess energy satisfies 5 > min (if 1 , If 5 )- 

The result just sketched shows that if 6 < min^ 1 , If 3 ) then for y € (—Zo, 0) the sheet is 
“fairly” confined. Though the rigorous result (Lemma 17.31) states this fact in an integral 
form, a heuristic version of its conclusion is that u x (— 1/2, y) ~ 1/2 and u x ( 1/2, y) ~ —1/2 
for y € (—Zo, 0). It follows that 

rl/2 rl/2 rl/2 

1 ~ u x (-l/2,y)-u x (l/2,y) < / \u x , x (x, y)\dx < / \u XjX +^ x /2\dx+ £, 2 x /2dx, 

J - 1/2 J- 1/2 J- 1/2 

so at least one of the terms on the right-hand side must be of order 1 (see Lemma 17.61 for 
the details). 

• The type I deformation considered in Section [6] involves a coarsening cascade of wrinkles. 
It is preferred over our other constructions in a certain (rather large) region of parameter 
space. Qualitatively, coarsening is favorable because coarser wrinkles have less bending 
energy. Lemma 17.71 makes this quantitative, by showing that the L 2 norm of the out-of¬ 
plane displacement £ cannot stay uniformly small. 
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• Continuing the discussion of coarsening: the boundary condition at the top makes £(-,0) 
small in L 2 , while it follows from the previous step that £(-,yo) is large for some yo € 
(—Zo, 0) (if the confinement is not released, and if the excess energy is small). In Lemma [775] 
we show that there is an energetic cost associated with changing the scale of the wrinkling; 
in fact, each doubling of ||£(-, y )\\ L 2 costs an energy of order Chy/rL. This lemma is, 
roughly speaking, a lower-bound analogue of the passage from (16.101) to (16.161) . Its proof 
builds on an argument from [25]. 

• The preceding bullets sketch the proof of the lower bound in the regime where self-similar 
coarsening is desirable. We also prove lower bounds associated with other regimes (for 
example when a type III deformation is best). If the bending term ||£, ra (-,y )||/,2 does 
not decrease much, we obtain the lower bound K 2 w^ 2 L by simply integrating the bending 
energy over the domain. Otherwise we can assume that for some yo € (—17/2,0) the 
bending term f £ 2 xx (x, yo) dx is much smaller than at the beginning (y = 0). In this case 
we consider y(x) := £(x, yo) — £(x, 0). We show that ||y,a:||L 2 is of order 1 while ||77,xx||z,2 is 
bounded by Cw^ 1 . These yield a lower bound for ||y||£ 2 , obtained as a consequence of the 
interpolation inequality \\r],x \\ 2 L 2 < C int (IMMI^xxlU 2 + ||r?||| 2 )- But ||y || L 2 is controlled 
by JJq £ 2 y , which in turn is controlled by the excess energy. We thus obtain another lower 
bound on the excess energy. 


7.2 The dichotomy 


We start with a lemma relating the excess energy to the out-of-plane displacement £ and the 
strain term u yiV + £ 2 y /2 — f yV . 

Lemma 7.2. Let (u, £) be any deformation satisfying (12.21) . and recall our definition of it excess 
energy: 

5 := E h (u,0 - B(f), 

where f is the minimizer of the bulk energy B (see Section W)- Under the assumptions of 
Theorem [7] the excess energy is nonnegative, and 


J \ u v,y + (iy/ 2 ~ f,y\ 2 dx < <*, // dy < 8 S(tL) 


n/2 


(7.2) 


JJ 2 \ Ux ,y + £,x£,?/| 2 + \ u x,x + £(L/2| 2 + /i 2 |V 2 £|“ dxdy < <5, (7-3) 

where fl/2 := (—1/2,1/2) x (— L/2, 0) is the upper half of the domain. 

Proof. We know that 

Eh{ u i £) = JJ \ u x,x + £/x/2| 2 + u x,y + Uy,x + £,x£,i/| 2 + \ u y,y + £,y/2| 2 + h 2 \ V^£| 2 dx dy 
+ t // Uy dx dy 


n 


> JJ \u y ,y+£,y/ 2 \ + TUydxdy = B(u y ) + JJ u y ^ y + ^ y /4dxdy. 

(7.4) 

Since / is the minimizer of B (in particular it is a critical point of B), using the previous relation 
we get 


1 


5 = E h (u, £) - B(f) > - ( D 2 B ■ (u y - f),u y - f) + JJ^u y! y£ 2 y + ^ y /4dxdy 

= Jlj u y-y ~ /. 2/) 2 + u v^l + CV 4dxd ^ = JJj u y-y ~ f<v + ^/ 2 ) 2 + f,y£l dxd v > °- 


(7.5) 
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We see that if E/(iq f) is close to B{f) then u y , y + t ^ y /2 is close to f. y and £ 2 is bounded. More 
quantitatively, we have 


/ \ u y,y + £«/ 2 “ f,y\ 2 d - Td l/ < £■ 

Jn 

T -^IL ildxiy -IL f ' A/2 - s ’ 

where we have used that f^ > 0 for y € (—L, 0) and f, y (y) > rL/4 for y € (—L/2, 0) (see 115.31) 1. 
Finally, we observe that in 117.51) we have used only one of the four non-negative terms in Eh, 
so (17.31) follows. □ 


We must consider the possibility that the sheet spreads laterally, releasing its confinement 
(at least in part). In the following lemma we estimate the energy needed to do that. 

Lemma 7.3. Under the assumptions of Lemma | 7. 2| there exist —1/2 < xq < —7/16 and 
—3/8 < x\ < —5/16, and a universal constant Cq > 0, such that either 

i) there exists l\ € (0, L/2) such that 

5 > Co 


and 


rx i 
J x 0 


u x (x,y) - u x (x, 0) dx 


< 1/512 


(7.6) 


for all y € (—Zi, 0); 

ii) or CHI) holds for all y € (—L/2,0) (we set l\ := L/2 in this case). 

Before proving the lemma let us sketch the main ingredients in its proof. First, we use the 
bounds from Lemma o to show that u y is in average close to /. In particular, f \u y (x,y) — 
f(y) | dy is small for generic x, and we can choose such generic xq and x\ near x = —1/2 such 
that the size of the interval (xq,xi) is of order 1. 

To get CHI) it is enough to estimate u X}V since f*' u x (x, y 0 ) -u x (x, 0) dx = // o ' /° Q u x>y dx dy. 
We write u x/y = (u XiV + u y , x + £, x €,y) ~ — £,x£, y - The first term on the right-hand side can 

be estimated using the excess energy 5 (see C3D). We integrate the second term in x and use 
that by the choice of xo,xi both u y (x o, •) and u y (x i, •) are close to /; it follows by the triangle 
inequality that u y (x o, •) is close to u y (x i, •). 

It remains to estimate £ t x€, y - First, we write £(x, y) = f £, y (x, -)+£(x, 0) and use Lemma[L2j 
Then using interpolation with £ jXX we obtain enough control over (, x . In the end we obtain an 

estimate for f^ 1 u x (x , y) — u x (x, 0) dx in terms of a function of <5max(Zi, if), which concludes 
the proof. 

In proving Lemma E3 and in many other places in this section, we will use the basic 
interpolation inequality 




(7.7) 


which holds for any interval I Cl and any (p € W 2,2 (I) 


with a fixed universal constant Ci n t- 
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Proof of Lemma [73| To start we need to choose xq and x\ generically. By (17.21) we have 

\u y ,y + i 2 y /2 - f,y\ 2 dx dy < 5, 


n/2 


£f y /2dxdy <46 (tL) \ 


so we can choose —1/2 < xo < —7/16 and —3/8 < x\ < —5/16 such that 

/ 0 rO 

\u y ,y(x i ,y) + (, 2 (x i ,y)/2-f,y(y)\ 2 dy<5 1 / ^ 2 (x i ,y)/2dy<5(TL)~ 1 , * = 0,1. 

-L J-L/2 

(7.8) 

Now suppose (USD is not true for all y € (—L/2,0); then there exists a smallest l\ E (0, L/2) 
such that 


rx i 

/ u x (x, Zi) Ux(X) 0) dx 

J XO 


> 1/512. 


By the boundary condition (12.21) we have u y (-,0) = 0. So (17.81) implies that for any y E (—Zi,0) 
and i = 0,1: 

I u y (xi,y) - f(y )| < f | u y , y (xi,y) - /, y (y)|dy 
Jy 

r0 rO 

< \uy t y(xi,y) + £, 2 y (xi,y)/2 - f, y (y)\dy + / ^(x*, y)/2 dy 
Jy Jy 

f \ /*o 

^ /2 v L l u ^( x * ,2/ ) + ^( a: * ,2/ )/ 2 -^( y )l 2dy J + J ^y( x ^y)/ 2d y 


Holder 

< 


117.81 

< 5 1 ' 2 i\ /2 + 5(tL)-\ 


We integrate the previous relation in y and use tL > 4 (see (13.11) 1 to obtain 

[ \uy(xi,y) - /(y)|dy < 8 1/2 l\ /2 + <5Zi, * = 0,1. (7.9) 

J-h 

We define J) 7 := [xq,xi] X [—Zi, 0]. Then 


rx i 

/ u x (x, -Zi) - ** x (x,0)dx 
J XO 


u X} y dx dy 


— II \ ^x,y T ** y,x T dx dy + 


Uy,x dx dy 


O' 


+ 


O' 


f.xf.y dx dy 


. (7.10) 


By (17.31) the first integral on the right-hand side satisfies 


Holder 




\ 1/2 


\'Ux,y T ^y,x T dx dy fZ H | ( jj T Uy }x T £,x£,j/| dx dyj ($ S ! Z^ 

(7.11) 


O' 


For the second integral we have 

r-0 


u y , x dx dy 


n 1 


< 


/ Uy(xi,y) - Uy(x 0 ,y)dy 

J-h 

[ \u y (xi,y) - f{y)\ dy+ f |%(x 0 ,y) - /(y)| dy < <5 1/2 Z^ /2 + <5Zi, (7.12) 
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where the last inequality follows from (17.91) . It remains to estimate the last term in (17.101) : 


f2' 


ix€,y dx d y 


Since we already have bound on £ ;y (see (17.21) 1 it remains to estimate £ jX . To do this, we first 
show a bound on £ by integrating then we use interpolation with £ )XX . For any y E [—h, 0] 
we have 


[ \£(x,y)\ 2 dx = f f £ i2y (x,t)dt + £(x,0) 

J xq Jx o Jy 

rx i rO 1 

<2 / / ^y(x,t)dt 

J xq J v 


dx 


'XQ Jy 

rx\ rO 


dx + 2 


rx i 

/ i£(®,< 


; dx 


< 


f*x~\_ rx\ 

2 \y\ / \£,y(x,t)\ 2 dt,dx + 2 |£(x, 0)| 2 dx < 51\{tL)~ 1 + u>o, 

Jy Jx 0 


where the last inequality follows from (12.21) and (17.21) . We use the interpolation inequality (17.71) 
for £ to get 


rx i 

I ICjxO^J 

j X 0 


dx < 


< 


rx\ \ / rx\ 

/ l£(-Ty)| 2 dx) f / |£, ra (x, y)| 2 c 

«/ iro / \J xq 

+ |xi -x 0 r 4 (y |£(x,y)| 2 dx^ 

(<5Zi(tL) -1 + Wq) f |£ )X x(x, y)| 2 dx + (<5 Zi(tL) -1 + wffi 

J x n 


After integration in y we obtain 


2 /*0 / rx\ \ 2 

2u^u„,l ^ 7 ' ' ' ''x,w)| 2 dx 


rO / rxi 

l£,x(x,y)| 2 dxdy J <h i |£ )X ( 

f2' / Jyo \Jx o 


</i((Hi(rL) +w 0 ) \ixx(x,y)\ dxdy+ h(5h{rL) + w Q ) 

JJn 1 

< h(8h(rL)- 1 + wl)5h~ 2 + h{5h(rL)- 1 + w 2 0 ) 2 , 
where the last inequality follows from (17.31) . This combined with (17.21) implies 


f}' 


£,x£,j/dxdy < | |C,ar| Iz, 2 (r2') I IC,x/1 IZy 2 (S7 / ) 

< ^(^2(rL)-i + <5h- 2 + /i(<JZi(rL)-i + (7.13) 

Combination of (|7.11D . (]7.12p . and (|7.13D with (|7.10p gives 

r x i 


rx\ 

/ w*(x, -/i) - 0) dx 

J XQ 


< 

r^j 


(S^l^ + d^l^ + dh] 

+ 5l\ /2 h- l / 2 {TL)~ 3/4 + ^/V/^/V 1 / 2 ^)- 172 + ^(rT)- 1 + tc 0 (5 1/2 /| /4 (rL) _1/2 . 

(7.14) 
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We want to estimate the four terms on the last line of (17.141) in terms of 8l\ and 8l \. By (17.11) 
we know 




<5 < C LB yiVrL log ( 


/ hi 


-i . .-1 


+ 1 I + WqtLI + l 


(7.15) 


\ WqVtL 

for any l € (0, L). It is easy to observe that for l > L the relation (17.151) holds as well since the 
right-hand side of (17.151) is then larger than CLBhyfrL log ^ ^ + 1^, which is in turn larger 

than 5 by (17.11) . We now estimate the last line of (|7.14l) by making three separate applications 
of (17451) : 

(1) Taking l := h _3 / 5 (rL) 1 / 10 in ()7.15jl and using that log(l + 1) < 2f 1 / 4 for t > 0 we obtain 


<5 < Clb Uh 3 ' 5 (J-'j V2 (rL) 2 / 5 + w 2 h 3 /\rL) 9 / 10 + h 3 /\rL)~ l ' 19 


(7.16) 


Since the first term in the last line of (|7.14|) can be written as 

8l\ /2 h~ l / 2 {rL )“ 3 / 4 = {81 3 ) 1 ' 9 (^-^(rL)- 9 / 10 ) 576 , 

we can control it by estimating <5 /i~ 3 / 5 (tL)“ 9 / 10 . Using h/wo < 1 (see (14.21) 1. rco < 1, and 
tL > 4 (see m), we obtain from (17.161) that 

5h~ 3 / 5 (rL)- 9 / 10 < C LB (2 (rL)" 5 / 10 + w 2 + (tL)" 1 ) < 1 , 

which implies 

<SZ} /2 /r 1/2 (rL)- 3 / 4 < ( 8l \) l/6 . (7.17) 

(The implicit constant in (17.171) involves a positive power of Clb • But we may (and shall) 
assume that Clb < 1; then (|7.171) holds with an implicit constant that is independent of 
Clb-) 

(2) For the next term we use (17.151) with l := h~ 4 . Using log(l + 1) < t 1 / 2 for t > 0 it follows 
from (17.151) that 

8 < Clb (+0 1 (rT) 1 / 4 + w\rLh + hj , 

so we have 


w 1 0 /2 8 3 / 4 l 1 1 /4 h~ 1 / 2 (rL)- 1 / 2 < {8k) 1 / 4 (8woh~ l {rL)- 1 ) 1 ' 2 

< C^i) 174 ((rL)" 3 / 4 + w 3 + w 0 (rL)- i y /2 < (8h^ 4 . (7.18) 


(3) The last two terms are fairly easy to estimate. First, we observe that 8(tL ) 1 < 1. Indeed, 
choosing l := 1 in (17.151) gives 




<5 < hVrL log 


/ h 




\wqVtL 


+ 11 + WqtL + 1 ] ^ 2 + WqtL , 


where we used log(l + t) < t and h < wo. Since tL > 4 and rco < c w (see (13.11) and (14. 2D ), 
we have 8(tL)~ 1 < 1 + ujq < 1. Using this estimate for 8, the last two terms are estimated 
as follows: 


«* /4 (rL)-‘ = (Shf‘‘(S(TLr l ) l, ‘(TL)- 3 /‘ < (Shf l ‘, 


where we have used (13.ID and (14. 2D . 
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Combining ([7.171) and (17.18j) with ()7.14j) and the last estimate we get that 


rx i 

/ u x (x , -ii) - ^(x, 0) dx 
J xo 


< (51 1) 1 / 2 + (.5Z 3 ) 1 / 2 + <5Zi + (<«?) 1/6 + (5Zi) 1/4 + (<5 Zi) 3/4 . 


By our choice of l\ we see that 


1/512 < 


rx i 

7 a?o 


u x (x,-l±) - u x (x,0)dx 


F(5 max(Zi, Z 3 )), 


where F(t ) = t 1 / 6 + f 1 / 4 + i 1 / 2 + t 3 / 4 + t. It follows that 

5 max(Zi, if) > Cq 


for a universal constant Cq > 0. □ 

Remark 7.4. By symmetry we can find, 5/16 < X 2 < 3/8 and 7/16 < X 3 < 1/2 for which either 
i) there exists I 2 E (0, L/2] such that 

5 max(Z2, Z|) > Cq 


and 


rx 3 

/ 

J X2 


u x (x,y) - u x (x, 0) dx 


< 1/512, 


holds for all y E (—Z 2 ,0) ; 

ii) or (17.191) holds for all y E (—L/2,0) (hre set I 2 := L/2 in this case). 
We use this remark together with Lemma 17.31 to get 


(7.19) 


Corollary 7.5. There exist —1/2 < xo < —7/16, —3/8 < xi < —5/16, 5/16 < X 2 < 3/8, 
7/16 < X 3 < 1/2, and a universal constant Co, such that either 

i) for any y E (-L/2,0) 



u x (x, 0) dx 


< 1/512, 


u x (x, 0) dx 


< 1/512, 


(7.20) 


ii) or there exists Iq E (0, L/2] such that 

5 max(lo, Iq) > Co 

and (|7.20l) holds for all y E (— Zo, 0). 

When the first alternative holds, we take the convention that Iq = L/2. 

Proof. We obtain this result by combining Lemma l7.3l with Remark [731 choosing Iq := min(Zi, fe). 

□ 

The next step in the proof of the lower bound is to show that f f 2 x (x, yo) dx and f \u XjX (x, yo)+ 
f 2 x (x,yo)/2\ dx can not be simultaneously small for any yo £ (—Zq, 0). 
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Lemma 7.6. For any yo E (— /o,0) we have either 

rx 3 

/ C 2 x (x,y 0 )dx > 1/4 


or 


rx3 

/ |«x,x(ab y 0 ) + ^(x, yo)/2| dx > 1/16. 

J XQ 

Proof. Let yo £ (—Zo>0) be fixed and xi E (xo,xi),x r E (x 2 ,xs). We have 

Us(av> 0) - u x {xi, 0) = , Z/o ) ~ u x {x u yo) + [{u x (x r , o) - u x (x r , yo)) - (u x (x h 0) - U x (xu l/o))] 


rx, 

Jxi 

rx, 

Jxi 


u x , x {x, Vo) dx + [(u x (x r , o) - u x (x r , yo)) - ( u x (xi , 0) - u x (x h y 0 ))\ 

1 1 f Xr 

u x , x (x,y Q ) + -i 2 x (x,y 0 )dx - - / £ 2 x (x, y 0 ) dx + [... ]. 

“ J Xl 


Therefore 


rx 3 2 r x r 

/ C%(x,y 0 )dx > - / ^ X (x,y 0 )dx 
J x o ^ J xi 

rx r J. 

> u x (xi,0) - u x (x r ,0) + / u XjX (x,y 0 ) +-£ 2 x (x,y 0 )dx+ [...] 

Jxi 
rx r 

> u x (xi, 0) - u x (x r ,0) - / 

J X 


u x , x (x,y 0 ) + ~^ x (x,y 0 ) 


dx + [... ]. 


We integrate the previous relation with respect to x/ and x r to obtain 


u x , x {x,y 0 ) + ~t.%(x,yo) 


^ rx 3 rx i rai3 r 

- £, 2 x ( x ,yo) dx > -r j- u x (xi,0)-u x (x r ,0)dx r dxi- 

^ J XQ J XQ J X2 J X>_ 

rx 1 rx 3 

- -f u x (x r , 0) — u x {x r , yo) dx r - -f Ux(xj, 0) - u x (x t , y 0 ) dxi 
J X 0 J X2 

i r 

i 2 x { x i 0) dx — 1/16 — / 

^ J X 1 J X< 


dx 


«x,x(®,yo) + ^, x (x,yo) 


dx, (7.21) 


where the last inequality follows from (|7.2U1) (averaging an integral adds a factor of reciprocal 


length to it, e.g. 


l l 


l l 


= t 4) and 


(which is used to replace u x (xi, 0) — 


512 |x 0 -xi| — 512 1/16 32 

u x (x r , 0) by f £ 2 x /2dx). We estimate the first term on the right-hand side by 

^ rx 2 ^ rx 2 r 5/16 

- / ^ 2 x (x, 0) dx = - / (2 cos(27rxu;(/ 1 )) 2 dx > 2 / cos 2 (27rxu? ( / 1 ) dx > 1/4, 

2 ./xi ’ 2 J Xl X—5/16 

where we used that f I cos 2 (9t)dt > |/|/4 provided that 0|/| > 2. (In the present setting 
|/| = 10/16 and 0 = 2-k/wq > 2-k/c w by 1)4.211 : since we always assume c w < 1, the condition 
holds with room to spare.) This concludes the proof since either we have 


rx 3 

J X Q 


u x , x (x,y 0 ) + -£ 2 x {x,y 0 ) 


dx > 1/16 


or else 


^ r x 3 ^ 1 f X3 

eAx ' m)i 


u x , x {x,y 0 ) + -£ 2 x(x,?/o) 


dx > 


□ 
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7.3 The energy required for coarsening 

When the type I deformation achieves the optimal scaling, we expect that the wrinkles must 
coarsen as one moves down the sheet. The following lemma justifies this, by giving a lower 
bound for the maximal amplitude of f. 

Lemma 7.7. There exists a universal constant C\ > 0 such that 

rX3 - 'h 2 l 0 


K := max 

ye(-lo,0) J XQ 


rx 3 

' m / f(x,y) 
>0) Jxn 


dx > C\ 


All. 


Proof. If 5 > lo/32 2 the assertion is almost trivial; indeed, since X 3 — xq > 1/2 and h < wq < c u 
(see we have 


K > 


£’ ?2(l ’ 0)d:c ~”»-'* 22 WTT - 


in this case. 

If on the other hand 5 < then we argue as follows. For any y € (—Zo, 0) Lemma 17.61 
implies 

rx 3 \ 2 rx 3 1 

j f 2 x(x,y)dxj +J \u x , x (x,y) + £ 2 x (x,y)/2\ dx > —. (7.22) 

Using the interpolation inequality (|7.7|) and \I\ 2 = \x% — xo\ 2 > 3/4 we obtain 




Hx(x,y)\ 2 dx 


' XQ 


< r z 

— °int 


px 3 \ I/ 2 / f x 3 \ I/ 2 Q fX 3 

/ \f,(x,y)\ 2 dx) I / \^ xx (x,y)\ 2 dx) +- \£(x,y)\ 

J XQ J \J XQ J ^ J XQ 

rx 3 

L / |£,3::e(x, y) \ 

J XQ 


dx 
2 dx + K 2 


We integrate (17.2211 in y and use the previous inequality to obtain 

'■0 / px 3 \ 2 /-0 rx 3 

—lo J xo 


/ 0 / rx 3 \ ^ /’O /- 3:3 

(/ ) dy+ / |u X)X (x, y 0 ) + £ )X (x, y 0 )/2| dx dy 

-Iq \Jx o / J—Iq J xq 

/ 0 rxs 

/ |C,a:x(x, y)| 2 dx dy + CIqK 2 

-lo j XQ 

/ r 0 f*XQ \ 

l o /2 ( / K,x(x,yo) +< 2 a: (x, 2 /o)/ 2| 2 dxdy j 

\J —lo J Xo J 

<f CK5h~ 2 + CIqK 2 + d 1 / 2 ^ 72 < CK5h~ 2 + CZ 0 A' 2 + l 0 / 32. 

using the assumption that 5 < Iq/{2>2) 2 . We have obtained 

Zo/32 < A5/W 2 + l 0 K 2 , 

which implies 


< CK 

+1 

m 


K = max 

ye(—/o, 0 ) J xq 


rx 3 

/ K(*,i/)| s 

J XQ 


dx > C\ 


h k 


A 1 


□ 
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The next lemma shows that coarsening costs energy: 


Lemma 7.8. Let 


. 1/2 


- 1/2 


a := I £ 2 (x,0)dx, b := max / £ 2 (x,y)dx, 

J- 1/2 S/e(-Jo,0) J-l/2 


and assume 


b > 4a, h 2 < Da (7.23) 

for some constant D > 0. Then there exists a constant C^iD) > 0, depending on D, such that 

tL rO ri/2 _ _ r ° r 1 / 2 

8 J—Iq J — 1/2 J-l 0 J-1/2 

r° r 1/2 

/ / | a x , x (x,y)+^ 2 c (x,y)/2| dx dy > C 2 {D)hVrL log(6/a), (7.24) 

4-/ n 4 - 1/2 


7° /' 1/2 

^(^yHxdy + h 2 / / £ 2 xa ,(x,y)dxdy 

4-Z n 4-1/2 


+ 


where b := b A (8Q n t) 

Proof. Let 

IV := Llog 2 (6/a)/2j > 1, 

where IV > 1 follows from (17.231) and the fact that 4a = 2u>Q7r” 2 < 2 c^7t~ 2 < (8Ci n i) _1 (the last 
inequality is a smallness condition on c w ). Then b > 2 2Jv a and for * = 0,..., N we can define 

.1/2 


Hi := max 


y € [—io, 0] : [ £ 2 (x,y) dx 

4-1/2 


= 2 2 *a 


(note that yo = 0). 

We will prove that for i = 0,..., N — 1 the left-hand side of (|7.24l) in (y ? ;+i, yf) x (—1/2,1/2) 
is of order at least hyfrL. Fix i £ {0,..., N — 1}. We see 


2 2i+2 a = 


rV2 .. 

/■1/2 

/ /•?/« \ 

II 

B 

+ 

N 

v/> 

/ 

/ GO^dy + £(®,yi) 

1 — 1/2 

/ —1/2 

\4 j/i+i / 


dx 


< 2 


< 


/■1/2 ryi r 1/2 

/ / C,y(^-2/)dy dx + 2 f 2 (x,yi) dx 

4-i/ 2 4j/ i+1 4 _i/2 


(7.25) 


2(yi-Vi+i) / y) dy dx + 

J — 1/2 J Ui -\-1 

which immediately implies 


2 • 2 2 *a, 


1/2 z-y; 

(yi-Vi+i) / ^ 2 y(x, y) dy dx 

J — 1/2 Jyi ~|-i 


> 2 2i a. 


(7.26) 


By Lemma 17111 for any y € (—Zo,0) either 

/■ 1/2 


'- 1/2 


|^x,s(ah y) + ^(z, y)/ 2| dx > 1/256 


or else 


ri/2 
' — 1/2 


f 2 (x,y)dx > 1/4. 
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In the latter case, for y E (jji+i , IJi) it follows from the interpolation inequality (17.71) that 

2 


1/16 < (j^J 2 x (x,y)dx\ 


< r * 2 

u int 


<2 Cf nt 


r 1/2 
7—1/2 
/• 1/2 
' — 1/2 


£ 2 (x,y)dx 

£ 2 (z,y) dxj 
1/2 


ri /2 

J —1/2 
- 1/2 _ 


\ 2 , 1/2 

y) dx j + J £ 2 (z,y) dx 


^LO^yjdx 1 + 2Cf nt b 2 


1/2 


< 2 C' 2 nt 2 2( * +1) a / £ 2 xa; (x, V) dx + 


2 C£ 


int 


- 1/2 64C?,’ 

and we can absorb the last term into the left-hand side to get 

- 1/2 


/ !/2 

£(L 0 uy)dx > !• 

1/2 


Using (17.231) we obtain 
- 1/2 
' — 1/2 


/ 1/2 ^ /* 1 /Z 

i“x,x(®, y) + £ %{x , y)/2| dx + h 2 / ^ xx (x,y) dx 
-1/2 7—1/2 

> min (l/256, h 2 / (a2 2 ( i+1 ))) > C(D) ^ 
We integrate this inequality in y over (yj+j, y,;) and use (|7.26l) to get 


2 2i a' 


tL f yi f 1 / 2 


ryi nr* ryi r 1/2 

/ / ^Ou^dxdy + h 2 / / £f xx (a;, y) dx dy 

«/7—1/2 7 y,’ _i_i J— 1/2 


Vi+l ^ 1/2 

/"m r 1 / 2 

+ 

^ Vi +1 / 1/2 


2 /i+l J — 1/2 
2 


rvi r 1/2 „ 9 T r o2i„ 

/ / i^x,x(x,y) + £ 2 x (x,y)/2| dxdy> —-b C(D) 

Jv<±. 1 7-1/2 8 y* - yi+i 


(y* - ym)// 2 


2 2 *a 

> C(D)hVrL 1 


where the last inequality follows from the AM-GM inequality. Since IV > log( 6 /a), summing 
the previous relation for i = 0,..., IV — 1 implies (17.241) . □ 


7.4 The lower bound half of Theorem |T| 


We are finally ready to prove the lower bound half of Theorem [TJ We assume throughout 
the following discussion that the assumptions of Theorem [T] are valid. As explained at the 
beginning of this section, we shall argue by contradiction; in particular, shall assume that (u,£) 
satisfies (EH). 

By Corollary 17.51 there exist —1/2 < xo < —7/16, —3/8 < x\ < —5/16, 5/16 < X 2 < 3/8, 
7/16 < X 3 < 1/2, such that either there exist Iq E (0, L/2) with the property that for any 
V € (—loi 0 ): 



u x (x,y) 

u x (x,y) 


— u x (x, 0 ) dx 

— u x (x. 0 ) dx 


< 1/512, 

< 1/512, 


(7.27) 
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and 


(7.28) 


5 max(Z 0 , 1$) > C 0 , 

or else (17.271) holds for all y € (—L/2,0) (in this case Iq = L/2). 
We use the notation from Lemma 17.81 i.e. 


a = 


f 1 / 2 „ 

/ C (®,o) 

J-1/2 


dx, 


6 = max 


" 1/2 


ys(-Zo,o) 7 _i/2 


£ 2 (x,y) dx A ( 8 C int 


.-I 


We distinguish between two cases: b > Aa (when the wrinkles are expected to coarsen) and 
b < 4a (when they are not). 

Case 1: b > 4 a. Using the same idea as in (|7.25jl we get that 

b<2 l 0 JJ £ 2 y dxdy + 2 a. 

Using b > 4a and a ~ iuq, the previous relation implies ff £ 2 y > uiq/q 1 . Now it follows from (17.21) 
that 


<5 > (tL) // £ 2 > wItLIq 1 . 


Next we want to show that 


5 > h\FrL log (w 0 2 [^= A l) + l) • 


(7.29) 


(7.30) 


We distinguish two subcases, when Mq/^Wq^/tL) < 1 and when HIq / {wqVtL) > 1. The defini¬ 
tion of the first subcase is equivalent to h\frL < WqtLIq 1 , which together with (17.291) implies 


5 > wItL1 0 1 > 


hVrL > hVrL log ( — + 1 ] > hVrL log ( w 0 2 ( 


hlo 


\wqVtL 


V u W^L 


A1 +1 , 


where the second-last inequality follows from the fact that the argument in the logarithm is less 
than 2 . 

Turning to the second subcase: we must show (17.301) when hlo / (w'qVtL) > 1. Let C\ be 
the constant that was chosen in Lemma rm We may suppose that 5 is small in the sense that 


5 < 2(CiTT 2 )h\frL log ( + 1^ , 

\w 2 VrL ) 


(7.31) 


since the opposite inequality implies (17.301) . By the initial conditions 

pl/2 


a = 


— 1/2 


£ 2 (x, 0 )d x = WqTt / 2 


in particular thanks to ico > h (see (14. 2D ), condition (I7.23D holds with D = 2n 2 . By Lemma 17.21 
25 is more than the left-hand side in (17.241) , and so using Lemma 17.81 and Lemma 17.71 we get 
that 

' 2 Citt 2 (h 2 l 0 


25 > C 2 hVrL\og(b/a) > C^hVrL log 


lOn 


A 1 


(7.32) 
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We now observe using (17.311) that 


2Ci7t 2 /i 2 / 0 hip 2C 1 T 2 hV^L 


w 2 5 


> 

w? 


o^xL 

hip _ 

i VtL i 0 g 


hip 


UqVtL 


+ 1 


SO 


C^VtT log ( 2C\ir 


,h 2 lp 




2 c I > C 2 hVrL log . 

WpdJ \ \ WnVrL 


hip 


- log log 


hip 


\VtL 


+ 1 


w; 


(7.33) 


Since log(t) — loglog(t + 1) > log(f + l )/2 for t > 1, combination of (17.321) and (17.331) implies 


5 > —j~hV tL min [ log 


h h \ .. (2Cm 2 

+ 1 ,4 log 


w. 


\VtL 


W n 


> hVrL log 


—2 


Wn 


(3H 


+ 1 


(7.34) 


where we used that hip/(wp VtL) > 1 and wp < c w (with c w sufficiently small). This completes 
the proof of (17.301) in the second subcase. 

We now easily conclude the validity of the lower bound half of Theorem [T] when b > 4a. In 
fact, if Ip < L/2 then adding (17.301) . (17.291) . and (17.281) gives 

6 > hVrL log (^Wp 2 == A 1 j + 1 j + w%t Lip 1 + min {Iq 1 , 3 } , 

while if Ip = L/2 we have 

5 > hVrL log 2 = A 1 ) + 1 


Case 2: b < 4a. Since b < 4a and 4a = 2wp7r 2 < 2 c^ 7 r 2 < ( 8 Ci n t) 1 for small enough c w , 
we immediately see that the value of b = max yg (_; o 0 ) Wi y) dx agrees with b. 

To prove the lower bound in this case it is sufficient to show that 


5 > min 


I 2 —2 

h w. 


L, 


mm 

l£(0,L) 


(h 


2 -2 

W. 


I + WqtLI 


-l 


+ min (^ 1 ,l 3 ))^ 


(7.35) 


Indeed, since t > log(l + 1) we have 

h 2 Wp 2 l > h VtL log ^ 9 ^__ + 1 ^ > h VtL log (^U-’p 2 A ^ + 1 

for any l > 0, and so the desired lower bound for 5 follows from (17.351) . 

The rest of this section is devoted to proving (17.351) . We note for future reference that the 
prescribed boundary conditions ( 12 . 21 ) satisfy 

!!£(■> 0 )IIl 2 (- 1 / 2 , 1 / 2 ) = W 0 7T /2 = a i 

IIC,x(')0)IIl2(_i/2,1/2) = 2) 

I |C,a:a:(') 0) 11L 2 (_i/2,1/2) = ^VT Wp . 
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If I2/)llia(_i/ 2 ,i/ 2 ) — 7r 2 ^o 2 /( 8C int) f° r al I V G (-^/2, 0) then our task is almost 
trivial, since 


171 

s > h z II £ 2 xx (x,y)dxdy>h 2 w 0 z L, 


(7.36) 


which implies (17.351) . 

Thus it suffices to consider the situation when ||£ iXa; (*,y )||^ 2 (_ 1 / 2 1 / 2 ) < 7t2 ' w o ~ 2 /( 8 Cf nt ) for 
some y E ( —L/2 ,0). Recalling (see Remark 1 7. II) that we may assume (it, £) are smooth, consider 
the first yo E (—L/2,0) (i.e. the one with smallest |yo|) such that 


rl/2 

I- 1/2 


£?**(*, yo) dx < tt 2 w 0 2 /( 8 C'f nt 


Note that by the definition of yo we have f £ 2 xx (x, y) dx > ir' 2 w 0 2 /(8C? nt ), for all y E (yo,0), 
hence 


m _ 

5 > h 2 


0 rl/2 


C 2 xx (x,y)dxdy>h 2 w 0 2 \ 


Let 

and observe that 
fl/2 


'yo -1-1/2 

y(x) := £(x,0) - £(x,y 0 ) for x E (-1/2,1/2) 


(7.37) 


\2 
IL 2 


— 1/2 


?y 2 dx = 


- 1/2 

-1/2 


l£(x, 0 ) — £(x, y 0 )| 2 dx 


< 


z-l/2 
' — 1/2 


I 2 / 0 I [ t 2 y (x,y)dy) dx < |i/o| ff t£ y dxdy. (7.38) 
J yo J »/»/ 


We now distinguish between two subcases: when |yo| > h and when |yo| < lo- Focusing 
initially on the subcase |yo| > ^o : our main task is to show that 


6>w%TL\ yo \ 1 . 


(7.39) 


Since yo 0 (—Zo, 0), the relation b < 4a does not imply J _^ 2 £ 2 (x, yo) dx < 4a. Therefore we 
must consider both of the following possibilities: 

(i) Suppose j/0 2 ^ 2 (x, yo) dx > 4a. Then the triangle inequality gives H^Hl 2 > 2yfa — y/a = 
fa = wo/{V2tt), and (17.381) gives 


^ydxdy < |y 0 |8(5/ (rL), 


w 2 0 /(2n 2 )<\\r,\\ 2 L2 < 
which implies (17.391) . 
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(ii) Suppose £ 2 ( x > 2/o) & x < 4a. Then we see using the interpolation inequality (17.711 

" 1/2 


HU-,2/o)|| i2( _ 1/2 , 1/2) = / C,x(*.yo)dx 


- 1/2 


/•V 2 ri/2 

/ C(x\y 0 )dx- / ^ x (a;,yo)da; 

J — 1/2 J- 1/2 


2 C 2 nt / 



1/2 


( f 1/2 

+ 2 C 2 nt 

/ i 


\J-l /2 

2 Ci n t 4fl 

TT 2 Wq 2 / 


i2 „2 


= 1/2 + 32C 2 nt a 2 < 1, 


where the last inequality holds for small enough a (hence small enough c w ). Therefore 
\\v,x\\l 2 > HC,x(-, 0)|| L 2 - ||^a;(-,2/o)||L 2 >V2-1, 

I I^xxIIl 2 — ||£,xx(‘) 0)1 Ii 2 — I |£,xx(": 2/0) I \l 2 w 0 • 

By interpolation (EZD 

(V2 — 1) < I|l/,x|l22 ^ C'int 111/11L 2 I |l?,xx 11L 2 T CWitlMIz, 2- 

The second term satisfies 

CintMh < C , i„ t (||e(-, 0)|| i 2 + ||C(-, 2 /o)||l 2 ) 2 < 1 6Cf nt a < SC^tT 2 < (y/2- lf/2 

for small enough c w . Thus it can be absorbed into the left-hand side of the previous 
relation. Combining these results with (17.381) . we get 


1 £ I2/0I ^JJ^ydxdyj w 0 2 , 

and so 

S>^ JJ^ 2 y > nftrLlyo]- 1 , 

confirming the validity of (17.391) . 

To conclude the treatment of the subcase |yo| > ^0 we must show that (17.351) holds. Com¬ 
bining (|7.37l) with (17.391) gives 

5 £ h 2 W Q 2 \y 0 \ + v^TL\y 0 \~ 1 . (7.40) 

If lo = L/2 then (17.351) follows from (17.401) (using that |yo| > Iq)- If 011 the other hand Iq < L/2 
then we have (using again that \yo\ > Iq) 

8 > h 2 WQ 2 \y 0 \ + w 2 TL\y 0 \-i + min(Z q Mo 3 ) 

> h 2 WQ 2 \yo\ + u!qtL|?/o | _1 + min(|y 0 | _1 , M -3 ), 

which implies (|7.35l) . The subcase |yo| > Zo is now complete. 

Turning now to the other subcase, when |yo| < lo : using Lemma 17.61 we see that for any 
y E (—Zo, 2 /o) either 

r 1 / 2 n 

1/4 < / (*.(x,y)dx 

4-1/2 
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or 


r 1/2 

1/256 < / \u XtX (x,y) + £, 2 x (x,y)/2\ 2 dx. 
J—1/2 


1 - 1/2 

If for a given y the first case is true, we use the interpolation inequality Cl 3 to get 

f 1/2 
/ — 1/2 


(7.41) 


/•V2 n 

1/4 < / C 2 a; (x,y)dx 
J — 1/2 


< Cint 


r 1/2 
1 — 1/2 


f(x,y)dx 


f 1/2 
-1/2 


£2 


\ ^ /- 1/2 

(x, y) dx j +C iri t J ?(x, y) dx. 


Since b < 4a ~ iuq and Cintda <1/8, the previous estimate implies 

r* 1/2 
1 — 1/2 

We combine this estimate with (17.411) to get 
ri/2 


/•i/2 

1^0/ £;L 0 uy)dx. 
J—1/2 


1 - 1/2 


\ux,xi x , 2/) + y)/2 | 2 + h 2 ^ xx (x, y) dx > min(l/256, Ch 2 w 0 2 ) > h 2 u ; 0 2 , 


/ 2Z0 cl/2 

/ K,*(®,y) + ? 2 c (x,y)/ 2| 2 + h 2 £ 2 ca; (x,y)dxdy > h 2 w 0 2 (l 0 - 12 /oI)- 

-Zn 7-1/2 


where we used that wq > h (see 621)). We integrate this in y over (—lo,yo) to get 
fyo ci/2 

/-z 0 7-1/2 
This combined with (17.371) implies 

5 > h 2 wfl 0 . (7.42) 

Since b < 4a, we have that f*y 2 £ 2 (x, yo) dx < 4a and arguing as we did earlier (in case (ii)) 
gives (rrm 

To conclude the treatment of the subcase |yo| < lo we must show that (17.351) holds. If 
Iq = L /2 then (17.421) gives 

^ > h 2 wfl 0 = h 2 wfL/2 , (7.43) 

which implies (|7.35l) . If on the other hand Iq < L/ 2, then we know 

<5 > h 2 WQ 2 l 0 + wlrL\y 0 \- 1 + min(Zo :1 , l ^ :3 ). 

Since |yo| < Iq, we have that uiQ-rLlyol ” 1 > WqtLIq 1 , this gives 

5 > h 2 wfl 0 + u^tLIq 1 + miri^ 1 , Zq 3 ), 

which implies (17.351) . This completes the treatment of the subcase |yo| < fo. The lower bound 
half of Theorem [T] has now been fully established. 
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